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I. 



INTRODUCTION 



A. BACKGROUND 

Jet engines undergoing maintenance are tested in jet 
engine test cells prior to their operational use. A typical 
test cell is shown in Figure 1 and can be viewed as a "U" 
shaped concrete structure of rectangular cross-section. The 
cell is designed to provide clean distortion free airflow 
to the engine inlet, which is necessary for accurate 
performance evaluation. In addition, the cell attenuates 
noise and provides safety for test equipment and personnel. 

Engine exhaust is directed into an augmenter tube, then 
deflected vertically to be expelled into the atmosphere. 

The augmenter tube acts in combination with engine exhaust 
momentum to form a jet pump. In reality, the engine exhaust 
jet entrains secondary air (air not ingested by the engine) 
which is then carried through the augmenter tube by exhaust 
jet momentum. Frictional losses associated with the mass 
of air moving through the test cell are thus overcome by the 
jet pump's energy addition. 

The quantity of secondary air entrained is crucial to 
proper engine testing. Too little entrainment allows exhaust 
gas recirculation to the engine inlet. Too much entrainment 
causes excessive static pressure gradients to exist between 
engine inlet and exhaust planes. Static pressure should be 
nearly identical in both planes for proper testing but a 
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small gradient must exist to prevent exhaust gas recircula- 
tion. Sufficient secondary air should be entrained to cool 
the engine exhaust to acceptable levels. Also, exhaust flow 
with very high kinetic energy can damage the test cell; 
secondary air entrainment reduces this energy. 

Augmentation ratio is defined as the ratio of entrained 
secondary air to the jet exhaust mass flow (primary flow) . 

In practice, engine/augmenter tube separation is adjusted to 
generate sufficient secondary flow to prevent exhaust gas 
recirculation at all power settings. Occasionally, the 
augmenter inlet diameter is varied to further control the 
entrainment process . 

Both analytical and empirical methods have been applied 
to the entrainment mixing problem. Empirical equations are 
usually applied to particular test cells, but operator 
experience coupled with trial and error are required to 
create the desired secondary flow. 

Larger and more expensive test cells are needed in the 
future to replace obsolete cells and to accommodate new 
engines with greatly increased mass flow rates. Present 
design tools do not adequately describe the flow field in 
a test cell and new design concepts cannot be evaluated until 
the cells are constructed. An accurate flow model would 
provide a needed design tool and would help prevent costly 
design errors. 
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B. TEST CELL POLLUTION 



Test cell installations are becoming increasingly sub- 
jected to localized state and community pollution control 
board regulations. The occasional prominent exhaust plume 
rising from test cells is a continual source of local com- 
plaint. Even when the plume is not visible, the cell emits 
pollutants similar to those of aircraft operating in the 
area. 

Fortunately, the stationary cell is amenable to exhaust 
treatment which has been shown to be effective in reducing 
particulate emissions to acceptable levels. The cost asso- 
ciated with the construction and operation of the necessary 
equipment is high and is a function of total mass flow 
(primary plus secondary) through the cell. Thus, it is 
necessary to minimize total mass flow to minimize abatement 
costs . 

Carbon particulates and unburned hydrocarbons are of 
immediate concern in the abatement of emissions from jet 
test cells. Also of concern are nitrogen oxides (NO ) , 
carbon monoxide (CO) and sulfur dioxide. The amount of each 
pollutant emitted is a function of engine type and mechanical 
condition. Power setting alters emissions with CO and unburned 
hydrocarbons predominate at idle. Particulates and nitrogen 
oxides are the main pollutants at high power settings . After- 
burner operation provides very high external exhaust tempera- 
tures and chemical reactions continue as the flow passes into 
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tiie augmenter tube. The secondary air mixing is believed 
to alter chemical reaction rates by effectively quenching 
portions of the exhaust jet. Present flow models do not 
describe the process completely and it has not been possi- 
ble until recently to accurately measure mass emission fac- 
tors in the afterburner core due to the high temperatures 
and the chemical kinetic nature of the process . 

To fully optimize the pollution abatement method it is 
necessary to understand how exhaust pollutants are formed, 
in what quantities and at what rates. The effect of engine 
power setting and secondary flow entrainment on emitted 
pollutants must also be analyzed. Pollutant measurement 
would be simpler and less costly if accurate species and 
velocity distribution predictions could be made at specified 
points in the test cell exhaust system. To accomplish these 
tasks, more complete knowledge of the air flow through test 
cells must be acquired. 

C. SCOPE OF STUDY 

Most analysis for predicting dilution air is very empiri- 
cal, usually requiring an empirical jet spread equation. 
Recent advances have been made toward solving the recircu- 
lating, turbulent fluid flow problem. One technique is due 
to Spalding et al. [Refs. 1 and 2], The governing two- 
dimensional elliptic partial differential equations are 
reduced to finite difference equations and solved by numeric 
methods . 
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Air flow through a test cell is known to be turbulent 
and does possess regions of recirculation. The purpose of 
this study was to apply the methods described by Spalding 
to the flow field encompassing the engine exhaust and 
augmenter tube of a typical Navy jet engine test cell. 

Flow visualization and static pressure recovery were 
primary objectives of this study. Parametric studies were 
made which involved the variation of augmenter geometry/ 
engine power setting, augmentation ratio and cell inlet 
flow conditions. 

The study was restricted to non-afterburning, subsonic 
engine exhaust conditions. 
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II. THEORETICAL APPROACH 



Two-dimensional fluid flow with regions of reverse-flow 
can be described mathematically by second-order elliptic 
partial differential equations. These equations can be 
reduced to "finite-difference" non-linear algebraic equations 
and solved simultaneously by an iterative procedure. One 
such procedure is detailed by Spalding and others in References 
1 and 2 and was the method used in this study. 

The conservation of mass, momentum and energy are the 
fundamental equations used in the problem formulation. The 
dependent variables chosen by Spalding et al. [1,2] are 
stream functions (t|0 and vorticity (w) to assure mass and 
momentum conservation. Temperature was used as the dependent 
variable for energy conservation since the flow was considered 
inert. Since the velocities under consideration were low, 
kinetic heating was ignored. Specific heat was considered to 
be constant. 

Although pressure and, to a large extent, velocity are 
removed from the formulation, they can easily be "backed-out" 
of the vorticity and stream functions solutions respectively. 

A turbulence model is needed in order to provide 
"effective" transport properties. The Jones-Launder 
turbulence model [2,3] is used to describe the variation 
of the effective viscosity, U e ff Two variables are used 
in this model: turbulence kinetic energy, K and turbulence 
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energy dissipation rate, e. The model is a time-averaged 
turbulence model and the two associated differential equations 
can be cast in elliptic form. Thus, the solution procedure 
for t(/, w and T also applies to K and e. The relationship 
between K and e is: 



^ef f 



C 

V 



P 




e 



where 



p = local density 

= empirically determined coefficient = 0.09 



These five equations together with the equation of state 
for a perfect gas describe the complete flow field. Choosing 
cylindrical coordinates (R,6,z) allows all five equations 
to be written in the following general form: 






aiJ b <f R ‘ 






+ R s 4, = 0 



where <p stands for ip, T, K, or e. The functions a^, 
b^, c^ and are given in Table I. 
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TABLE I 



COEFFICIENT VALUES 



<J> 



a b 



c 



s 



0 ) 

R 



eff 



S (to/R) 



* 



K 



e 



0 



1 



1 



P R 



eff 

a K 



eff 

a 



1 



1 



1 



to 

R 



- (G - pe) 



-( 



c^Ge 

K" 



P ^ 



2 

-) 



T 



1 




1 



0 



where 



y 

p f f 

a v 5 r taken to be 1.0 

K j K,eff 



r K,eff 



the exchange coefficient of the kinetic 
energy of the fluctuating motion, K. 
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a 



E 



v eff 

r e,eff 



taken to be 1.3 



r e , ef f 



the exchange coefficient of the turbulence 
energy dissipation rate, e. 



o 



H 



M eff 

F h,eff 



taken to be 1.0 



?h eff the exchange coefficient of heat 



G 



= y 



ef f 



\ 2 [ 

i. 



3V 

+ 



9V R,2 



,V 

v R ' 



+ 




+ 




c^ and C 2 are empirically determined coefficients. 
( c ^ = 1.45, c 2 = 2.0) 



Neglecting tangential momentum, the vorticity source term 



S (o,/R) 15 



S (u/R) - - R 



L 92 



( 



v 2 + V 2 
z V R . 



9R 



- ( 

9R v 



V 2 + V 2 
2 R x i£ 
2 ’ 9z 



+ negligible terms 
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The first term of the general equation is the convective 
term while the second and third terms are the diffusion terms. 
The last term is the source term which includes all terms not 
strictly convective or diffusive. 

Once the equations are solved and a converged solution 
is acquired, the static pressure distribution in the flow 
field can be calculated by integrating the momentum equation 
in the axial or radial direction. The resultant equation for 
the gradient in the axial (z) direction is 



The analogous form in the radial (R) direction (neglecting 
0-momentum) is 




1 3 

R 3r 







where 
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Pressure is then obtained by the following approximate 



integration formula!: 



B 

/ 

A 




1 t( i£) 

2 U 3z ; 



(*£) 

l 3z ; 



B 



(Z B - 



Z A> 
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III. MODEL DISCUSSION 



A. COORDINATE SELECTION 

Although the test cell cross-section is basically rectan- 
gular, the engine exit and augmenter tube can be accurately 
described using two-dimensional axi-symmetrical flow. Of 
primary interest is the mixing of the engine exhaust (primary 
flow) and ambient air (secondary flow) . This process pri- 
marily occurs between the two circular "tubes" and within 
the augmenter tube. Reference 4 suggests that test cell 
wall proximity will influence primary/secondary flow mixing. 
However, the large dimensions of the cell compared with the 
engine and augmenter diameters should allow the use of 
cylindrical coordinates with little error. 

B. GEOMETRY 

Figure 1 is a profile view of a test cell located at 
NARF Alameda, used to test TF-41 engines. Figure 2 is the 
geometry established for an engine test while Figure 3 is the 
geometry used in this study. The actual dimensions are 
easily varied in the computer program allowing parametric 
studies to be made to assess the effects of geometry change 
on the flow field. The field is assumed axi-symmetric although 
in practice the engine/augmenter combination is located 
nearer the test cell floor than the walls or ceiling. The 
cell ceiling is modeled as cylindrical but in actuality is 
quite flat. 
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To minimize the number of grid points required in the 
axial (z) direction, the cell "inlet" was located a few feet 
upwind of the engine exhaust. 

The straight augmenter duct "lip" used in the computer 
model was a necessary simplification to the beveled lip 
since grid points in the coordinate system lie along mutually 
perpendicular grid lines. Augmenter tubes in Navy test cells 
are by no means uniform in design and inlet shapes vary 
greatly. It is believed that a straight lip design does 
not greatly affect the flow analysis. 

Augmenter tubes used to test afterburner operation 
usually contain quench systems and flow baffles. These were 
not considered in this analysis. 

C. GRID SPACING 

Various grid sizes were used to appraise the effect of 
grid size on the predicted flow field. The largest number 
of grid points was described by a 60x50 (zxR) matrix and the 
smallest by a 20x25 matrix. Regardless of the size of the 
grid, every boundary must be defined by a grid line. In 
general, grid line separation should be smaller as one 
approaches a wall since gradients in properties are large 
in this region. 

Spalding [1,2] suggests that in the region of solid 
boundaries each succeeding grid spacing not differ from the 
last by more than a factor of 1.5. 
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D. BOUNDARY CONDITIONS 



All dependent variables must be specified in some manner 
on the boundaries. The boundaries of this model include: 

(1) the engine exhaust duct, (2) the cell inlet, (3) the 
solid boundaries, (4) the augmenter outlet, and (5) the 
geometric centerline. 

1 . Engine Exhaust Duct 

a. Stream Function/Velocity 

Plug flow was assumed at the engine exhaust 
plane. The axial velocity distribution is uniform (Figure 
4) to produce a specified mass flow and no radial component 
of velocity exists. 

b. Vorticity 

Since no velocity gradients exist, vorticity 
is identically zero. 

c. Mass Fraction/Temperature 

A non-dimensional mass fraction/temperature was 
assigned a value of 1.0 at the engine exit. This corresponds 
to the actual exhaust temperature existing at the specified 
power setting. 

d. Turbulence Kinetic Energy 

Kinetic energy of turbulent motion (K) is defined 
by 



K 





+ V 



) 



where 



v 



( ) 



are the instantaneous fluctuations 
in free stream velocity. 



24 



For isotropic turbulence the expression for K 



3 T 

can be written K = ^ v . 



Intensity of turbulence is defined as 



I = 



J 



1 , 2 . 2 . 2 . 

3 V R * v 6 + v z > 

V 



Assuming isotropic turbulence 



I = 






V 



where V is the mean flow velocity. 

Spalding [1,2] presents an empirical relation 
for K in uniform flow: 



K = 0.004 (V) 2 



Setting 



3 T — 2 

j v = .004 (V) z 

and substituting into the intensity relation we find 

7^008 v 2 

I = 4 : = 0.052 

V 
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Schlichting [5] experimentally determined 



intensity values of between 0.05 and 0.08 in low subsonic 
wind tunnel flow. Thus, it appears Spalding's relation is 
adequate for estimating K in the engine exit and was utilized 
as the boundary condition. 

e. Turbulence Energy Dissipation Rate 

K and e are related through the dissipation 
length scale Z [2] by 

k 3/2 
e = 

Z 

In this program Z was taken to be: 



2 . Inlet to Cell 

a. Stream Function/Velocity 

A linear variation of inlet axial velocity was 
chosen which varied from a positive value at the engine 
wall to zero at the cell ceiling to produce a specified 
secondary mass flow (Figure 4) . No specific information was 
available to specify this flow distribution accurately. 
Lacking more specific information, the radial components of 
velocity (V ) were fixed at zero. 

x\ 



Z = 0.022 x (Engine Radius) 



b. Vorticity 



Vorticity was computed using 



w 

R 
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c. Mass-Fraction/Temperature 



Non-dimensional temperature and mass-fraction 
were assigned a value of 0.0 at the cell inlet. This 
corresponded to an antual temperature of 520 °R. 

d. Turbulent Kinetic Energy (K) 

K was specified at each grid point using the 

expression 



K = 0.004 V 2 

_2 

Where V is the local flow velocity, 

e. Turbulence Dissipation 
K ^/l 

0.022 x (Cell Radius - Engine Radius) 

3 . Solid Boundaries 

There are eight solid boundaries described in the 
present model. Each is numbered in Figure 5 for reference. 
Two are engine surfaces, two are cell walls and four describe 
the augmenter tube. 

a. Velocity /Stream Function 

Along the engine surfaces (1 and 2, Figure 5) 
stream function is constant and determined from the specified 
mass flow through the engine. The remaining solid boundaries 



e = 
l = 
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have the stream function fixed at a value which is the sum 
of the engine mass flow and cell inlet mass flow. All 
velocities at a solid boundary are zero, 
b. Vorticity 

Vorticity near a solid boundary was computed in 
one of two ways depending on the magnitude of the local grid 
cell Reynolds number. 

(1) Laminar Sublayer Case. If the local grid 
cell Reynolds number (R = Aip/yR) is found to be less than 
290 [1,2], a linear velocity gradient is assumed. Using the 
additional assumption that the velocity gradient in the 
direction perpendicular to the wall is much greater than the 
gradient parallel to the wall, the vorticity is calculated 
using the vorticity definition and discarding the parallel 
term. For example, given a wall parallel to the horizontal 
(z) axis with flow in the positive z-direction, the vorticity 
would be found by: 



0) 



3Z 




3V 


AV 


V 


z 


z 


- _£ 


3R 


AR 


AR 



3V R 

where -r — is considered small. V is the velocity in the 

<3 Z ir 

z-direction at node P (Figure 6) . Since velocity at the wall 
is zero, it is easily seen that V p is twice the average 
velocity between the wall and the nearest grid point off the 
wall. This is convenient since the average velocity is 
related to the stream function by 
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V 

z 

avg 



V 

P _ 1 Aif> 

2 pR AR 



and the local Reynolds number is 



Re 



y 



AijJ 

R 

P 



Therefore, the vorticity expression is given by 



to 





2 Re 
p (AR) 



which is easily calculated. 

(2) Turbulent Case. When the local grid cell 
Reynolds number at the near wall node exceeds 290 the flow is 
considered turbulent. Vorticity is then defined as 



to = 



T shear 

'Vail 



where r , is the wall shear stress and y . . is the 
shear p wall 

viscosity on the wall. T s hear 9^- ven by the relation: 



shear 



Cf V 2 
-y p V 



The coefficient (-^-) was calculated using 



.0284 

“ _ :T5 

Re 



where Re = 



Ajj) 

y K 
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(3) Special Points. The preceding discussion 



of vorticity calculation near solid boundaries assumed a 
well-developed boundary layer. There are locations in the 
field where the fluid may meet a parallel solid surface and 
the vorticity calculation becomes suspect. Such a point is 
the intersection of surface 5 and surface 6 of Figure 5 
(point C) . There are two problems associated with this point. 
If the flow is predominantly in the positive z-direction, a 
fully developed boundary layer will not exist on surface 5 at 
point C. If the flow is in the negative R-direction near 
surface 6, the identical problem arises. Also, each surface 
(5 and 6) will have a unique value of vorticity at point C 
and it is necessary to keep them properly identified in the 
computer program (see sketch below) . 



• • 0 

surface 5 

* r — 

c 



surface 

6 



D 



i 



R 




z 



No completely satisfactory numerical method 
exists in this scheme to account for the singularity. Initial 
computer runs were made allowing the vorticity at point C 
to be calculated as if a fully developed boundary layer 
existed on both surfaces (5 and 6) . Converged solutions 
resulted with realistic flow visualizations around this point 
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as long as a recirculation region did not develop on the 
augmenter duct lip (surface 6) . 

When flow developed in the positive R-direc- 
tion near surface 6 with negative z-direction flow on surface 
5, the solution would not converge and often diverged (see 
sketch below) . 



Surface 

6 




References 1 and 2 state that special 
attention must be paid to points which define enlargements 
in pipe flow (Points A and D of Figure 5) . In order to 
obtain realistic jet spread, stream function must be fixed 
on the node immediately downwind of these points. It became 
evident that the flow would usually be in the negative z- 
direction along surface 5, so ip was fixed at the node circled 
in the sketch above. This solved the divergence problem 
occurring under certain flow conditions, and appeared not to 
alter the solution appreciably. Point D is unusual in that 
it defines the end of an "infinitely thin" wall. It also 
defines an "enlargement" and for this reason, the immediate 
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downward ijj value was fixed at the value assiqned to 
point D. 



surface 

6 




D 

• • 



F 

> 

dominant flow 

direction 



For purposes of calculating the vorticity at node F (see 
above sketch) , the vorticity at point D was taken to be the 
average of the vorti cities computed on surface 6 and 7 shown 
in the above sketch. 

c. Turbulence Kinetic Energy 

Turbulent flow produces steep gradients of all 
transport properties near walls . To reduce the required 
number of grid points and to reduce computation time a 
"slip value" of turbulence kinetic energy (K) is calculated 
on the wall. The "slip value" is fixed to produce an arith- 
tnetic-mean gradient between the wall and adjacent node, which 
equals the gradient of the Couette flow solution midway 
between the nodes. A complete discussion of this method can 
be found in References 1 and 6 . 
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d. Turbulence Dissipation 

The effective boundary for turbulence dissipa- 
tion (e) is taken to be the n grid points adjacent to the 
wall. Wolfshtein [6] presents the value at the near-wall 
node as : 



e 



P 



(^4 (K p ) 3/2 /AR 

y 



where 




near-wall value 

empirical constant = 0.09 

empirical constant = 0.22 

near-wall turbulence kinetic 
energy value 



AR = wall to near-wall node 
distance 



e. Temperature/Mass Fraction 

All solid boundaries were assumed adiabatic 
and, of course, impermeable. Therefore, no gradients in 
temperature or mass fraction existed between the walls and 
near-wall nodes . 

4 . Augmenter Exit 

Axial gradients of all dependent variables were 
taken as zero. 

5 . Geometric Centerline 

Radial gradients of w/R, K, e, T and mass fraction 
are taken as zero. Also the centerline is the stream function 
(<Jj = 0) reference line. 
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E. INITIAL CONDITIONS 

The flow field must have initial values assigned to all 
variables. Experience is required to fill the field with 
compatible variables to prevent divergence in the first few 
iterations . 

It was found that a reasonably accurate stream function 
distribution was essential to prevent divergence. In par- 
ticular, the predominant flow direction must be established 
and large discontinuities in stream function should not 
exist. The flows in the augmenter tube and the free stream 
area between the engine exhaust and augmenter tube were 
initialized as well-developed pipe flow. This allowed com- 
patible vorticity values to be calculated in these regions . 

Density, mass fraction, viscosity, and turbulence 
kinetic energy were specified using reference values taken 
from either the cell inlet or engine exit. Turbulence 
dissipation was calculated using the relationship discussed 
under boundary conditions . 

F. SOLUTION PROCEDURES 

To counter computational instabilities in the first 400 
iterations, stream function variation between iterations was 
under-relaxed by a factor of 0.3. All other variables were 
under-relaxed by 0.5. The relaxation parameters were 
increased to 0.6 and 0.8, respectively, after 400 iterations. 

Various flows, geometries, and grids were investigated 
successfully using the above relaxation parameters. They 
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precluded divergence in the most unstable flow condition 
handled: high engine mass flow with low augmentation ratio 

(200 lbm/sec, 0.25, respectively). However, it should be 
noted that initial conditions must be well established for 
these relaxation parameters to be effective. 
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IV. MODEL DEVELOPMENT 



The geometry shown in Figure 3a was chosen to represent 
a test cell section containing an engine exhaust nozzle and 
augmenter tube. Approximate dimensions were taken from a 
TF-41 turbofan engine undergoing testing at NARF Alameda. 

A. GRID SPACING 

A grid composed of twenty lines in the z-direction and 
twenty-five lines in the R-direction (20x25) was chosen to 
initially describe the flow field. (Intersections of 
mutually perpendicular grid lines define grid points ♦ ) 

This number of grid lines allowed the use of small grid 
spacing near solid walls but produced relatively large 
spacing away from walls. 

Computation time is directly proportional to the number 
of grid points in the field. A small number of grid lines 
was chosen to minimize computer time during program debugging. 
Three hundred iterations typically required 7.5 minutes 
(excluding compile time) . It should be noted that vorticity 
source terms were ignored during this initial period. 

Vorticity source terms and pressure recovery calculations 
were added only after the program was successfully iterating 
to a converged solution. The number of grid lines was 
increased to 40x40 (zxR) to allow for pressure gradient 
calculations without the use of wall or near-wall variable 
values. Otherwise, large errors accumulated in the pressure 
integration process. 
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Eleven hundred iterations typically required 109 minutes 
of computation time (excluding compile time) . This time 
includes one pressure calculation along eight lines in the 
flow field and those computations performed in the flow 
initialization process. However, the total time for the 
latter calculations amounts to less than one minute. 

B. CONVERGENCE CRITERIA 

Convergence rate was influenced by flow conditions, grid 
spacing and relaxation parameters . Two criteria were used 
to appraise the convergence level of the solution. 

First, the percentage change in hj, K, e and mass 
fraction (at all grid points in the field) must be arbi- 
trarily small between successive iterations. This was done 
computationally by storing the largest change occurring in 
the field for each variable and printing that value after 
each iteration. 

The second convergence criterion proved to be the most 
conclusive. That is, when pressure calculations along each 
grid line remained unchanged (with increased iterations) 
convergence was assumed. 

C. PLOTTING PACKAGE 

Solution values of stream function and temperature were 
transferred to computer cards at the conclusion of each run. 
These cards were then used as input to a second program 
which plotted streamline contours and lines of constant 
temperature on separate graphs. A CALCOMP Model 765 plotter 
was utilized for this purpose. 
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V. PARAMETRIC STUDY 



Geometry and inlet conditions were varied to appraise 

model utility and sensitivity to flow conditions. TF-41 

engine test data were used to establish engine exhaust static 

temperature and velocity. Static pressure at the engine 

2 

exit was fixed at atmospheric pressure (2116.0 lb f /ft ) for 
this study. Engine exit velocity was found by computing p 
from the equation of state for a perfect gas and using the 
one-dimensional mass conservation equation. 



P 



P 

RT 



and 



V 



engine 



m 

engine 
p A 

K engine 



where m 
A 



total mass flow through engine 
(including bypass air) 

area of exhaust plane 



Engine exit static temperature was computed using the 
turbine outlet temperature and the bypass ratio given in 
Reference 7. A summary of these parameters is shown in 
Table II. 
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TABLE II 



TF-41 TEST CONDITIONS 



Power 

Setting 


m (lb m /sec) 


T exit <° R) 


V exit (ft/sec) 


Exit 

Mach 

Number 


Idle 


100 


800 


594 


0.43 


80% 


150 


900 






90% 


200 


927 


1376 


0.92 


94% 


238 


955 






97% 


249 


970 






100% 


263 


1001 







A specific flow condition was chosen to provide a solution 
to which all succeeding solutions could be compared. That 
condition was: idle engine power setting with an augmentation 

ratio (AR) of 0.5. The geometry was fixed as shown in 
Figure 3. 

The "reference" flow condition was chosen as it provided 
low subsonic flow conditions at the engine exhaust plane 
(M = 0.43) and the augmentation ratio was typical for actual 
test conditions. 
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Pressure calculations were made along the dashed lines 
shown in Figure 7. The lines show only approximate locations 
since grid spacing was often altered when geometry was 
varied. 

Eight computer runs were made to examine the usefulness 
of the model. Each run was allowed approximately 2.5 hours 
(1400 iterations) to assure solution convergence although 
most runs achieved practical convergence prior to this 
time. A summary of the study follows in Table III. 



TABLE III 
FLOW CONDITIONS 



1. Reference flow condition was established 

as discussed above: Idle Power, AR = 0.5, 

Figure 3a geometry. 

2. Two runs were made varying the AR: Idle 

Power, AR = 0.75 and 1.0, Figure 3a 
geometry . 

3 . Two runs were made to change the length 

of the augmenter duct lip: Idle Power, 

AR = 0.5, Figure 3b, geometries. 

4 . One run varied the engine to augmenter 

separation distance: Idle Power, AR = 0.5, 

Figure 3c, geometry. 

5. One run varied the augmenter radius: 

Idle Power, AR = 0.5, Figure 3d, geometry. 

6 . One run changed the power setting of the 

engine: 90% power, AR = 1.0, Figure 3a 

geometry . 
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VI. RESULTS AND DISCUSSION 



A. PLOW CONDITION 1 (Reference Condition) 

Streamlines and static temperature contours are shown 
in Figures 8a and 8b, respectively. Pressure gradients along 
axial and radial grid lines (Figure 7) are graphically 
displayed in Figures 8c through 8f. Pressure contours along 
axial and radial grid lines are displayed in Figures 8g 
through 8i. Figure 8j presents the variation in augment er 
pressure rise as a function of iteration number. 

1 . Streamlines/Temperature Contours 

Axisymmetric flow solutions are obtained with this method. 
Therefore, boundary conditions on stream function are 
determined from 

pV R = " (1/R) 9l|,/8Z and pV Z = (1/R) 3l|,/8R 
For example, the mass flow rate of a TF-41 engine is 100 lb m /sec 

at idle power setting. Then, for uniform inlet velocity and 

density 

Ait = i|/ - ip = ip = T7— = = 15.9 lb /sec-radian 

y r r o 2ir 2ir m 

as is shown in Figure 8a. 

Two distinct regions of recirculation are apparent 
in Figure 8a. One prominent region exists above the augmenter 
tube bounded on two sides by the test cell ceiling and wall. 
Points of stagnation are shown located on the exterior 
augmenter tube surface and cell wall. 
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The second prominent region of recirculation is 



located downwind of the augmenter inlet lip. Close exam- 
ination of the streamline having its origin at the lip 
reveals axial flow for a short distance followed by rapid 
directional change. Axial flow could have been the result 
of fixing the stream function at the first downwind corner 
as discussed in Section III .D. 3 .b. ( 3) . However, an additional 
computer run was made with the same flow conditions except 
that the stream function at the downwind node was allowed to 
be computed normally. The solution for stream function was 
essentially identical to the first run. In any case, the 
lip has only small effects on the shape and size of the 
recirculation region . The recirculation region is formed 
at the low thrust, low augmentation ratio condition pri- 
marily by the radial component of the augmentation air as 
it flows into the augmenter. 

The region of flow reattachment is well defined on 
the augmenter interior surface and is located approximately 
two augmenter diameters downstream of the inlet. 

Temperature contours shown in Figure 8b are an 
indication of the amount of mixing occurring between 
engine flow and secondary air. From the non-dimensional 
mass fraction and temperature definitions, a thoroughly 
mixed flow would have a static temperature of 707 °R. 

Under these flow conditions, a well-mixed flow occurs 
approximately four augmenter diameters downstream 
of the augmenter inlet. The temperature contours 
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clearly show the effect of secondary air entrainment and 
recirculation within the augmenter tube on the temperature 
field. 

2 . Pressure Recovery 

a. Axial ( z-direction) 

Pressure gradients along lines A-A through E-E 
(Figure 7) are shown in Figures 8c, d and e. Lines labeled 
A-A and B-B are axial lines having their origin at the 
engine exit plane and terminating one node upstream of the 
augmenter exit. Lines C-C, D-D and E-E extend between the 
third node in the axial direction and the node immediately 
upstream of the augmenter exit. 

Large negative pressure gradients exist along all 
lines from the engine exit plane to just inside the augmenter 
inlet plane. Lines A-A and B-B exhibit the largest gradients 
in this region and examination of the solution reveals the 
exhaust jet rapidly departing from the uniform inlet velocity 
(plug flow) condition. Secondary air is being entrained in 
this area with line B-B being nearest the turbulent mixing 
zone of the exhaust jet. This may account for the largest 
pressure gradients existing along line B-B. 

Pressure gradients along axial line F-F (Figure 7) 
are very small and are not shown. Static pressures along 
axial lines A-A through E-E are shown in Figures 8h and i. 
These pressure profiles are all initialized at 2116.0 
lb f /f t 2 . 
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b. Radial (R-direction) 

Pressure gradients and static pressures along 
G-F are shown in Figures 8f and 8g, respectively. The engine 
exhaust jet turbulent mixing zone is clearly indicated by 
fluctuating pressure gradients (Figure 8f) . Static pressure 
is seen to increase toward the engine exhaust jet centerline 
which is due primarily to the secondary air momentum in 
the negative R-direction . 

No appreciable pressure gradients existed along 
radial line H-H which is located downstream of the flow 
reattachment point. 

c. Total Axial Static Pressure Rise 

Total rise (AP) was defined as the static pressure 
difference between end-points of the lines discussed above. 

AP is plotted against number of iterations in Figure 8 j . The 
upwind end pressure for all lines is reference pressure but 
each line is seen to exhibit a unique value for AP at 1400 
iterations. No pressure gradients exist in the radial 
direction downwind of the reattachment point. Therefore, 

AP values should be identical in magnitude for lines A-A 
through E-E when the solution is converged. 

Close examination of the pressure gradients 
(Figures 8c through 8e) reveals line E-E to have small, "well 
behaved" gradients compared with lines A-A through D-D. 

Figure 8j indicates that the AP calculation for line E-E 
"converged" relatively early (1000 iterations) . For these 
reasons, pressure rise was taken to be that computed for 
E-E (47 lb f /f t 2 ) . 
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The differences in the pressure rise calculated 
along axial grid lines can be due to several things . The 
number of iterations has a dominant effect on the calculated 
pressure rise. The dependent variables (w, \p, K, e, T) vary 
only slightly with increased iteration time. However, the 
pressure calculation is apparently very sensitive to small 
changes in the stream function and density. This is probably 
because is a function of second derivatives of \p. In 

addition, pressure rise is calculated by assuming linear 
variations in gradient between grid points and pressure has 
been held fixed for purposes of calculating the density 
(i.e., p = p ( T) ) . Improved accuracy might be obtained if 
pressure was calculated periodically in the solution in 
order that pressure effects on density may be taken into 
account. It was found that finer grid spacing did not 
appreciably affect the pressure calculation. Another possible 
error source may enter through the use of upwind differences 
and/or the inadequacy of the Jones-Launder model for the flow 
geometry under investigation. However, checks on the 
solution indicated that mass was conserved at each axial 
location with less than 14 percent error. 

B. FLOW CONDITION 2 

Augmentation ratio (AR) was varied with idle power setting 
to assess the effect on total static pressure rise through the 
system. Solution convergence and pressure gradients were 
similar to Flow Condition 1 and axial line E-E was again 
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chosen to determine the value of AP. The solutions along 
the other axial lines had approximately the same slope on an 
AR vs. AP plot but were different in magnitude as discussed 
above. 

Figure 9a shows decreasing AP with increased AR. This 
result qualitatively correlates with Bailey [8,9] who experi- 
mentally varied exit pressure on a small jet pump and deter- 
mined secondary mass entrainment. As the augmentation air 
is increased for a fixed engine flow there is less momentum 
differential between the two streams. This apparently causes 
the decreased pressure rise with increased AR. In actual 
test cell operation there is a unique AR for each engine 
flow rate (assuming the augmenter is fixed in design and 
position) . The actual value is determined from the pressure 
rise in the augmenter tube being equal to the pressure losses 
through the rest of the test cell. If pressure losses 
through the inlet and exit stacks (with their associated 
acoustic and pollution abatement equipment) and through the 
horizontal portion of the cell can be estimated as a function 
of flow velocity, then the AR can be found in an iterative 
manner by using the AP vs. AR plot determined with the model 
for a particular augmenter design and position. 

Figures 9b and 9c show the effect of AR on the stream 
function and temperature distributions, respectively. As 
the AR is increased the axial velocity component of the 
secondary air dominates over the radial component. The latter 
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component was observed above to be the major cause of the 
recirculation region. Thus, as the AR is increased, the 
recirculation region within the augmenter is seen to decrease 
in size and mixing occurs more rapidly. 

C. FLOW CONDITION 3 

The augmenter duct lip length was varied by 0.63 ft 
(2.1 ft to 2.73 ft augmenter inlet radius). Stream function 
and temperature fields are shown in Figures 10a through lOd. 

No appreciable change in calculated pressure rise was observed. 
This result apparently occurs because at low augmentation 
ratios the size of the inlet lip has very little effect on 
the size of the recirculation region within the augmenter. 

A greater effect of duct lip length on pressure rise could 
be expected at higher augmentation ratios. 

D. FLOW CONDITION 4 

Engine to augmenter inlet separation was increased to 
2.5 ft (from 1.5 ft) . Stream function and temperature 
contours are included as Figures 11a and lib. No appreciable 
change in pressure rise was noted along axial grid line E-E 
when compared with Flow Condition 1. 

The stream function solution indicated that no appreciable 
change occurred within the augmenter when the distance was 
increased. Again, a more significant effect could be expected 
for higher augmentation ratios. Moving the augmenter closer 
to the engine should have a major effect on the recirculation 
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region and the calculated pressure rise. As discussed above, 
the inlet velocity profile for the augmentation air was 
assumed to increase linearly from the cell ceiling to the 
engine at a location just upstream of the engine exhaust. 

This profile would change as the augmenter position is changed. 
Since test cells are designed to provide uniform flow at the 
engine inlet, future investigations with the model should be 
made with the secondary air inlet specified as uniform at the 
engine inlet plane. In this case, the velocity profile near 
the engine exhaust plane would be free to develop in a 
realistic manner. 

E. FLOW CONDITION 5 

Augmenter radius was increased by 0.5 ft to 3.5 ft while 

keeping the duct inlet diameter fixed. Figures 12a and 

12b present the temperature and stream function contours. 

Pressure rise (AP) along E-E was diminished to 37 lb^/ft 

2 

from the reference condition value of 47 lb^/ft . Although 
only one run was made, this pressure loss is considerable an 
indicates AP decreases with increasing augmenter diameter 
(given fixed inlet flow conditions) . The increased size of 
the recirculation region (which consists primarily of 
augmentation air) may contribute to this pressure loss. 

F. FLOW CONDITION 6 

Engine power setting was increased to 90% and the A R 
was set at 1.0 which is typical for an actual engine test. 
Figure 13a shows the recirculation region smaller with flow 
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reattachment farther upstream. The temperature contours 
(Figure 13b) show engine exhaust and secondary air being 
thoroughly mixed farther down the augmenter tube than 
occurred at idle flow. 

Pressure rise versus number of iterations is shown in 
Figure 13c. The numerical values of the dependent variables 
did not change appreciably after 1400 iterations. However, 
AP's differed in value considerably along lines A-A through 
E-E, so no definitive pressure rise could be derived from the 
results. The solution procedure (which assumes incompressible 
flow and uses the associated elliptic equations) used in 
this study is questionable when applied at the higher subsonic 
Mach numbers . 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The model provides axisymmetric flow visualizations 
in turbojet test cells and augmenter tubes for low subsonic 
flow conditions. These can be used to identify regions of 
recirculation and to assess the amount of mixing occurring 
between engine exhaust gases and secondary air. Optimum 
locations for pollution sampling equipment can be selected 
by examining the numerical solutions. 

Static pressure rise within the augmenter tube can be 
calculated for a given flow condition and augmenter design 
from converged solutions, provided that regions of "steep" 
gradients in variables are avoided. Total pressure rise 
through the system can be found and factors which influence 
this rise (augmenter design and location, engine flow rate, 
test cell size, and etc.) can be studied. 

High subsonic flows invalidate the use of partial 
differential equations of the elliptic type; therefore, 
the model provides realistic solutions only at low subsonic 
Mach numbers . 

B . RECOMMENDATIONS 

1. Additional computer runs should be made to firmly 
establish the effects of augmenter and cell geometry on flow 
mixing and pressure rise (AP) . For example, the effects of 
moving the augmenter inlet closer to the engine exit plane 
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and moving the test cell wall forward to the augmenter inlet 
plane should be considered. In addition, the program could 
also be easily modified to extend the augmenter duct lip in 
the opposite direction. 

2. The cell "inlet" position should be moved upstream 
to the engine inlet plane by the addition of a few grid 
lines in the z-direction. Plug flow could be assumed at 
the "new" cell inlet and the resultant flow field would show 
the actual velocity distribution at the "cell inlet" used 

in this study. In addition, this would permit calculation 
of exhaust gas recirculation. 

3. Static pressure and velocity measurements should be 
taken under various test conditions at NARF Alameda. The 
size and shape of the recirculation region located above 
the augmenter should be qualitatively verified by photo- 
graphing a tufted grid placed between the test cell augmenter 
tube and ceiling. These data would help in validation/ 
improvement of the computer model. 
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FIGURE 1. TEST CELL GEOMETRY 
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FIGURE 2. TF-41 TEST CONFIGURATION 
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FIGURE 3a. COMPUTER MODEL GEOMETRY: 
FLOW CONDITION 1 
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FIGURE 3b. COMPUTER MODEL GEOMETRY: FLOW CONDITION 3 
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FIGURE 3c. COMPUTER MODEL GEOMETRY: FLOW CONDITION 4 
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FIGURE 3d. COMPUTER MODEL GEOMETRY: FLOW CONDITION 5 
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FIGURE 6. LAMINAR SUBLAYER VELOCITY PROFILE 
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FIGURE 7. LOCATION OF GRID LINES FOR PRESSURE CALCULATIONS 
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FIGURE 8a. STREAMLINE CONTOURS: FLOW CONDITION 1 
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FIGURE 8b. TEMPERATURE CONTOURS: FLOW CONDITION 1 
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FIGURE 8c. AXIAL PRESSURE GRADIENTS: FLOW CONDITION 1 

(LINES A-A AND B-B) 
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FIGURE 8d. AXIAL PRESSURE GRADIENTS: FLOW CONDITION 1 

(LINES C-C AND D-D) 
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FIGURE 8e. AXIAL PRESSURE GRADIENTS: FLOW CONDITION 1 

(LINE E-E) 
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FIGURE 8f. RADIAL PRESSURE GRADIENTS: 
FLOW CONDITION 1 
(LINE G-H) 



radial 




FIGURE 8g. RADIAL STATIC PRESSURE VARIATIONS: 
FLOW CONDITION 1 
(LINE G-H) 
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FIGURE 8h . AXIAL STATIC PRESSURE VARIATIONS: 
FLOW CONDITION 1 
(LINES A-A AND B~B) 
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FIGURE 8i . AXIAL STATIC PRESSURE VARIATIONS: FLOW CONDITION 1 

(LINES C-C THROUGH E-E) 
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FIGURE 8j . AP RISE VERSUS NUMBER OF ITERATIONS: 
FLOW CONDITION 1 



— i 

1500 



66 




FIGURE 9a. AP RISE VERSUS AUGMENTATION RATIO: 
FLOW CONDITION 1 
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FIGURE 9b. STREAMLINE CONTOURS: FLOW CONDITION 2 (AR=1.0) 
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FIGURE 9c. TEMPERATURE CONTOURS: FLOW CONDITION 2 (AR = 1.0) 
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FIGURE 10a. STREAMLINE CONTOURS: FLOW CONDITION 3 

(SMALL AUGMENTER INLET RADIUS) 
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FIGURE 10b. TEMPERATURE CONTOURS: FLOW CONDITION 3 

(SMALL AUGMENTER INLET RADIUS) 
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FIGURE 10c. STREAMLINE CONTOURS: FLOW CONDITION 3 

(LARGE AUGMENTER INLET RADIUS) 
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FIGURE lOd. TEMPERATURE CONTOURS: FLOW CONDITION 3 

(LARGE AUGMENTER INLET RADIUS) 
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FIGURE 11a. STREAMLINE CONTOURS: FLOW CONDITION 4 
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FIGURE lib. TEMPERATURE CONTOURS: FLOW CONDITION 4 
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FIGURE 12a. STREAMLINE CONTOURS: FLOW CONDITION 5 
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FIGURE 12b. TEMPERATURE CONTOURS: FLOW CONDITION 5 
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FIGURE 13a. STREAMLINE CONTOURS: FLOW CONDITION 6 




79 



FIGURE 13b. TEMPERATURE CONTOURS: FLOW CONDITION 6 
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